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Abstract. 

Consider two identical atoms in a spherical harmonic oscillator interacting with a 
zero-range interaction which is tuned to produce an s-wave zero-energy bound state. 
The quantum spectrum of the system is known to be exactly solvable. We note that 
the same partial wave quantum spectrum is obtained by the one-dimensional scale- 
invariant inverse square potential. Long known as the Calogero-Sutherland-Moser 
(CSM) model, it leads to Fractional Exclusion Statistics (FES) of Haldanc and Wu. 
The statistical parameter is deduced from the analytically calculated second virial 
coefficient. When FES is applied to a Fermi gas at unitarity, it gives good agreement 
with experimental data without the use of any free parameter. 
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1. Introduction 

Experimental investigations on ultra-cold gas of fermionic atoms near Feshbach 
resonance, in recent years, have opened new avenues to address and understand the 
problems in strongly correlated fermionic systems pQ. Two identical fermionic atoms 
trapped in different hyperfine states may still interact in the relative s-state. Low 
energy properties of such a gas at low density are determined by the scattering length 
a, the number density n and the temperature T. The effective attraction between the 
atoms near a Feshbach resonance may be increased continuously by varying an applied 
magnetic field to reduce the Zeeman splitting between the states occupied by the atoms 
and the resonance. The scattering length a goes from a small negative to a positive 
value. The unitary limit is achieved when \a\ is infinite at the transition when a changes 
its sign resulting in a zero-energy two-body bound state. This defines the unitary limit 
and the behaviour is expected to be universal (scale independent) in this limit [2]. 

Recently, Liu et al [3] have calculated the virial expansion coefficients of the equation 
of state of a strongly correlated trapped Fermi gas on either side of the unitary limit, 
extending the work of Ho and Mueller [4j. The latter had earlier developed a virial 
expansion up to the second virial coefficient to study the universal behaviour of a 
homogeneous gas at unitarity. The central point of these investigations, for our purpose, 
is the fact that the second and third virial coefficients, when plotted as a function of the 
interaction parameter (scattering length in this case), become temperature independent 
at unitarity, testifying to the universal nature of the quantum gas at this limit. More 
over, the two-body bound state spectrum (in the s-state) for harmonic confinement 
at unitarity, shown in Fig. ([I]), exhibits the striking property of an overall shift in the 
energy levels due to the interaction. This is a hallmark of the inverse square interaction, 
which, in one-dimension, also leads to fractional exclusion statistics(FES ) as defined by 
Haldane [5|. In FES, the occupancy factor n»(T) of a single-particle state with energy 
Cj at temperature T is given by [6] 

rii = (wi + g)~ 1 , (1) 

where the distribution function w satisfies the nonlinear relation 

wfil + Wi) 1 -* = ia)/3] . (2) 

In the above, g > is the (temperature-independent) statistical parameter, (3 = 
l/(fcsT), and \x the chemical potential. The parameter g is based on the rate at which 
the number of available states in a system of fixed size decreases as more and more 
particles are added to it. As such, g assumes values and 1 for bosons and fermions 
respectively, because the addition of one particle reduces the number of available states 
hyg ®. 

As may be deduced from Eqs. (0Q EJ) the occupancy factor rij at T = for an 
ideal FES gas is specially simple, and is given by n« = 1/g up to fcj < hp, and zero 
otherwise, where hp is the shifted Fermi wave number. The relationship between hp 
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Figure 1. The s-wave spectrum of few lowest states of two identical fermionic atoms 
are shown: (a) Spectrum of non-interacting fermions. (b) spectrum of interacting 
fermions at unitarity. Note the spectrum is simply shifted down by one unit of energy. 
The energies are given in units of huj, and the spacing between two-adjacent s-wave 
states is 2~hu. 



and the fermionic kp is obtained by noting that the particle number N is 

1 r^F Air 1 ~ 

N = — Ank 2 dk = —-kl . (3) 
g Jo 3 g F K ' 



But for fermions N = ^-k F . For a fixed N, we thus get k F = g 1 ^k F . A similar 

r an iHdqI T7TV.^s rrac miroc Tr. — 

g 5 2M 



1 A fc 5 ~ 

calculation for the energy E of an ideal FES gas gives E = --f^ff ■ Eliminating kp, we 



then get 

E_ = 3^4 

N s 10M ' v ; 

where £ = g 2 ^ 3 . 

In an earlier investigation [8] , the energy per particle and the chemical potential at 
finite temperature of the quantum Fermi gas at unitarity were calculated by mapping 
the interacting fermionic system to a system of non-interacting quasi-particles obeying 
FES. In a subsequent paper [9], some properties of few-body systems were calculated 
in the same scheme. The statistical parameter g was determined phenomenologically 
from the energy per particle of a unitary Fermi gas at T = 0, given by Eq. (j3J). A 
strongly interacting Fermi gas at unitarity has no length scale other than the inverse of 
the Fermi momentum. Consequently, its potential energy has the same /cp-dependence 
as that of the kinetic energy. Since the potential is attractive, the parameter £ is less 
than unity. The parameter £ = 0.44 is close to its experimental value [10] and therefore 
the statistical parameter is g = £ 3//2 = 0.29. Since the value of the statistical parameter 
is dependent only on the nature of the interaction, it is fixed once and for all at all 
temperatures. Using this value of g and the distribution for FES particles given by 
Wu [6], the average energy as a function of temperature was calculated. The agreement 
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with the Quantum Monte-Carlo calculation (QMC) [TTI [T2"] for a homogeneous gas was 
found to be satisfactory. The idea was then extended to harmonically trapped gases 
where it was found to agree not only with earlier calculations [13] . but also with the 
available experimental data [H] (See Fig. 2, which is taken from [8]). The approach in 
our earlier work was phenomenological, involving just one scale independent parameter 
g. In FES, it has been shown [15] that g is determined by the high-temperature limit 
of the second virial coefficient. Our purpose in this paper is to determine the statistical 
parameter g from the second virial coefficient. The second virial coefficient is obtained 
from the inverse square potential in terms of g using a semiclassical procedure, which 
is known to reproduce the exact quantum result [16J. Liu et. al have also obtained it 
directly from the quantum spectrum given in Fig. ([I]). Equating these two second virial 
coefficients, we find that g = 1 — l/y/2 ~ 0.29. With this result, our FES calculations 
done earlier [8] require no free parameters any more. 




Figure 2. Energy per particle as a function of temperature. At left we give the results 
for a homogeneous gas (solid line). Our results are compared with the MC calculations 
of ref. [12] (solid squares) and ref. [IT] (solid triangles). On the right the results for a 
harmonically confined system is shown (solid line) . The dashed line corresponds to the 
calculations presented in ref. [13] and the experimental data are from ref. [14]. See [8] 
for more details. 

The relative s-wave two-body spectrum of Fig. ([I]) was calculated for two interacting 
atoms in a spherical harmonic oscillator by Busch et. al [18J. The interaction is of zero 
range, and its scattering length is tuned to infinity. The spectrum is universal, and is 
practically unaltered for interaction potentials whose range is much shorter than the 
oscillator length. In section 2, we briefly recapitulate the essentials of this spectrum. 
In section 3, we show that this spectrum is also reproduced by an equivalent one- 
dimensional system, namely the Calogero-Sutherland model (CSM) [17], with an inverse 
square interaction. This lends credence to the connection of the unitary spectrum to 
FES, since it is well known that CSM provides an exact realisation of FES. Section 4 
contains the main theme of the paper. Here g is obtained from the high-temperature 
second virial coefficient of the gas, as outlined earlier in this section. In FES, once 
the temperature parameter g is fixed, the properties of the gas are determined at all 
temperatures. Our calculations show that not only is the cold-atom two-body problem 
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at T = mimicked by the two-body Calogero model, but the same applies for the 
many-body case at finite T. 



2. Two-body spectrum at unitarity 

Consider two identical fermionic atoms in different spin states interacting with a zero- 
range potential. When the strength of the interaction is tuned to produce a zero- 
energy bound state, the scattering length a = ±oo. When these two atoms are trapped 
in a three-dimensional spherical harmonic oscillator, the full spectrum is analytically 
calculated and is well-known [T8l [19| [20] . The s-wave relative spectra, after subtracting 
the Centre-of-Mass (CM) energy, for the lowest few states are shown in Fig.flT]). It is 
important to note that the only effect of the interaction is to produce a constant unit 
downward shift in the energy levels, which are labelled in units of the oscillator spacing 
hu. As pointed out earlier, this behaviour is also a characteristic of CSM [17J with an 
inverse square interaction. The un-normalised ground state eigenfunctions in relative 
coordinates, uo(r) = r^o(r), for the two cases are given by 

Uo(r) = r exp(— r 2 /2), non-interacting, 

Uo(r) = exp(— r 2 /2), interacting, r > . (5) 

where the relative distance r is dimensionless, expressed in units of the oscillator length 
L = yj (h/mu), and m = M/2, M being the mass of the atom. The tower of states built 
on the ground state are the nodal excitations. Note that in the interacting case the wave 
function actually corresponds to the irregular solution of the non-interacting system that 
is normally excluded as an eigenstate. However, this is a valid solution at unitarity due 
to the presence of the singular interaction at the origin [20]. The spectra shown above 
remains almost unchanged even when the interaction range is finite, provided the latter 
is much smaller than the oscillator length [20J. 

In the next section we illustrate a one dimensional template wherein such a 
spectrum is realised. 



3. s-wave spectrum from Calogero model 

We now demonstrate that the two-body s-wave spectra shown above may be generated 
by the one dimensional Calogero Hamiltonian [IT]. We start with the two-body 
Hamiltonian 

where the interaction strength is controlled by g > 0. We may transform to CM co- 
ordinates 

P=(pi+p 2 ), X = (xi + x 2 )/2; 

and relative co-ordinates 



V = (Pi - P2)/2, x = (xi- x 2 ). 
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The CM Hamiltonian is that of a single oscillator and does not play any further role. 
Consider the relative Hamiltonian given by 

The particles cannot cross in one dimension and we may restrict the solutions to the 
region < x < oo (this mimics the range of r in three dimensions). Because of the 
singular nature of the interaction, the solutions go to zero at the coincident point. We 
define r = x/L, and express the energy in units of Hcu. The physically acceptable ground 
state solution in the interval < r < oo and its energy are again given by 

ipa(r) = r 9 exp(— r 2 /2), 

E =(g + 1/2) . (8) 

The full spectrum of states is easily found, with the energy eigenvalues and the 
corresponding eigenstates given by 

E n =(2n + £ + l/2), n = 0,l,2,3,.. (9) 

Mr) = r 9 exp(-r 2 /2)L^ 1/2 (r 2 ) . (10) 

Furthermore, following Calogero, the physically acceptable solutions may be extended 
to the whole range — oo < x < oo by imposing the condition 

ip(-x) = ±ij;(x) . (11) 

From the above symmetry/antisymmetry condition, we note that g — 1 corresponds to 
noninteracting "fermions" , and g = to "bosons" . This interpretation in one dimension 
should not be taken literally, since the permutation of the particles by crossing is not 
allowed. Comparing Eqs.flH]) and (El), the spectra (a) and (b) in Fig.l, we see that (a) 
corresponds to g = 1, and (b) to g = 0. It is important to note that the interaction 
vanishes at both g = 1 and g = 0, but the tower of energy states are not identical. In 
the first case, g = 1, the states correspond to the non-interacting system shown in the 
s-wave spectra of the three dimensional system in Fig. ([I]). The interaction is attractive 
for g < 1 and repulsive for g > 1. The maximum attraction is precisely at g = 1/2. 
In the second case, g = 0, is approached from the attractive side and the spectrum is 
identical to the interacting case shown in Fig.([T]). 

It has been shown by several authors [21] that the interacting particles of CSM may 
be regarded as non-interacting quasi-particles obeying FES with g as the statistical 
parameter as defined by Haldane [5]. In this sense the spectra shown in Fig.(TjQ) are 
remarkable. In the light of CSM we may physically interpret the s-wave spectrum 
obtained using the pseudo-potential as similar to the phenomenon in which the 
interaction is statistical in nature which produces the effect of turning fermions into 
bosons. While this analogy is indicative of the nature of the interaction as statistical, 
the actual value of the statistical parameter for the which the spectra agree, namely 
g = in one dimension, cannot be interpreted literally in three dimensions. 
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4. The universal second virial coefficient 

The grand partition function of a system may be expanded as a series in fugacity 
parameter, z = exp(— /3/i), at high temperatures (or low densities). The second virial 
coefficient a 2 = —b 2 , where 62 is the so-called cluster integral [22]. The second virial 
coefficient, to a large extent, determines the thermodynamic properties of a dilute 
interacting gas. If in particular the interaction is statistical in the sense defined by 
Haldane, the second virial coefficient is related to the statistics parameter g in the 
high temperature limit and plays an important role in determining systems which obey 
FES [15]. A classic example is the Calogero model where a gas with inverse-square 
pairwise interaction can be regarded as an ideal gas obeying FES |21j . 

The contribution to the interacting part of the second cluster integral may be 
expressed in terms of the spectra of interacting and non-interacting two particle systems 
and is given by 

A6 2 = b 2 -b° 2 = S n (exp -(J3E n ) - exp(-/3£°)), (12) 

where E n (E®) corresponds to the relative energy of the interacting (non-interacting) 
system. Substituting for the spectra as shown in Fig. 1, we then get 

Ab2 = (^tm. _ fHLW (13) 

\l-exp[-2y) l-exp(-2y)/ 
_ exp(-y/2) 



1 + exp(-y) ' 

where y = Koj(3. Taking the high-temperature limit y — > 0, we get 

+ (15) 

Although calculated in harmonic confinement, the temperature- independent value of 
1/2 is universal, and is also valid in a homogeneous gas. This is in agreement with the 
result of Beth and Uhlenbeck |23j , and Ho and Mueller [I] who obtained the universal 
value of 1/2 at resonance for a homogeneous gas. 

In the above, we have not included spin factors. If the spin factors are included 
in the one-body canonical partition function we have Ab 2 = 1/4 in agreement with the 
pseudo-potential calculation of Liu et al |3j. The overall factor of 1/2 due to spin will 
be omitted also in the subsequent semiclassical calculation to be consistent. We exploit 
this universal value of A62 at the s-wave resonance to provide a microscopic explanation 
of the origin of FES in cold fermionic atoms. The CSM type inverse-square interaction 
is scale invariant and gives rise to a temperature-independent universal second virial 
coefficient in any dimension. We therefore use the semiclassical approach to provide 
a link between Ab 2 and strength of inverse square interaction. In the partial wave 
decomposition we have 
1=00 

A6 2 = £(2Z + 1) A$\ (16) 
1=0 
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where the contribution due to the interacting part, that is Ab%\ can be written in the 
semiclassical WKB approximation as 



Ah 



(0 



dr exp 



h 2 i(i + i) 

Mr 2 



[exp(-/5V(r))-l], 



(17) 



where A = y 1 zttTi 2 f3 / M is the thermal wave length, V(r) is the two-body potential. 
Furthermore summing over all the partial waves, treating / as a continuous variable, we 
get 

2tt 



Ab, 



x r 2 dr [ eX p(-^(r))-l], 

which is indeed the correct semiclassical expression [22] . In general, the classical Ab^ 
in Eq. fll7l) even at resonance depends on temperature as well as the parameters of 
the potential. The lowest order WKB approximation is poor at resonance, and the 
universality is lost. The only exception to this rule is the scale invariant inverse square 
potential. More over, when the Langer modification [TB] of replacing 1(1 + 1) by (Z + l/2) 2 
is implemented, the WKB approximation reproduces the quantum results exactly. As 
seen from the one dimensional example, the s-wave asymptotic wave function is exactly 
reproduced at resonance when the scattering length a — > oo. Note however that the 
inverse square potential is only applicable in the / = partial wave at resonance. The 
inverse square potential is a long-range potential, and one may ask why its contribution 
to the second virial coefficient from the I > channels are not included. To answer this, 
recall that we are calculating the interaction part of the second virial coefficient. Due 
to the FR in the I = partial wave, the interaction is dominant only in this channel, 
and the higher partial waves contribute negligibly. In Busch et al.'s paper [18], the two- 
body spectrum is is obtained from a pseudopotential that acts only in the s-state. This 
spectrum is also used by Liu et. al [$]. At unitarity, there is no length scale due to the 
interaction, and we take the effective interaction to be inverse square, only applicable 
in the s-state. Away from the unitary point, the fractional exclusion statistics (FES) is 
not applicable. 

We therefore assume that the two body s-wave potential in relative coordinates is 
given by 



V (r) 



n 2 g(g 



M 

Substituting this in Eq.([[ 

get 

1 

A 



>. (19) 
r~ 

setting 1 = 0, and implementing the Langer correction, we 



Ah 



(o) 

2 



/■oo 




/ dr exp 
/o 


\-e " 2 1 

H AMr 2 



[eM-PV (r)) - 1] 



V2 l 2 

Equating this to the universal value of A62 
Eq.ffToT). we have 

1 

2 



(20) 
(21) 



Ab, 



1 r 1 r K 



obtained from the s-wave spectrum in 



(22) 



Cold atoms at unitarity and inverse square interaction 



9 



We have two solutions corresponding to the g — 1 — 1/ y/2 ~ 0.29 and g = 1/ \[2 =~ 0.71. 
The solution g = 0.29 is appropriate for an attractive interaction in the fermionic basis. 
Note that this result would not change by taking the spin into account, as both sides of 
Eq. (1221) would be multiplied by 1/2. An identical result is obtained when instead of a 
homogeneous interacting gas, we put the particles in an oscillator trap. Now Eq. (l20l is 
given by 

[exp(-(3V (r)) - exp(-^Wr 2 )] , (23) 







/ dr exp 
/o 


\-p e 1 



V (r) = tM^ + ^^L_ll. (24) 



where 

-Mu 2 r 2 + — 
4 Mr 2 

The interaction parameter^ may still be interpreted as the statistical parameter of FES 
since the mapping between the inverse-square interaction in one dimension (CSM) and 
FES is exact. The result of the integration in Eq. (l23~]) is 

1 



exp(-huj(3J(g - 1/2) 2 ) - exp(-hu(3/2) 



(25) 



In the limit of hu/3 — > the result is the same as in Eq. (!2TT) . It may be noted that the 
harmonic oscillator merely acts as a regulator in the high temperature limit, and yields 
the homogeneous gas result. 

The relation between the second virial coefficient and the statistics parameter g for 
a homogeneous gas in FES [15] is given by 

\-9 = 2 d/ %, (26) 

where d is the dimension of the space that is relevant. For s-wave contribution alone, 
we choose effectively d — 1. We therefore have 

l-g = V2Ab 2 , (27) 

where the non-interacting limit has g — 1 for fermionic atoms. This reproduces the 
solution given above for the interaction parameter g in CSM thus establishing the 
connection with statistics parameter of FES. 

The ground state energy of a gas of FES particles is given by Eq. (plj. Even though 
this relation is for a three dimensional gas, we use g obtained from the one-dimensional 
relation (1271) to determine £. This is because the deviation of the parameter £ from the 
noninteracting value of unity is due to the attractive potential energy. This potential 
energy arises from the pair-wise interaction that (as already noted) acts only in the 
I = state. So far as the potential energy is concerned, only a single partial wave is 
relevant, and the system is effectively one-dimensional. 

We have seen that the value of g obtained above from the high-temperature regime 
is compatible with the experimental results at all temperatures, thus establishing the 
connection between an ideal gas obeying FES and a dilute gas of strongly interacting 
fermionic atoms at unitarity. The earlier phenomenological analysis of the average 
energy and the chemical potential at finite temperature using the distribution function 
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for FES particles is also therefore justified. As further evidence, we may add the 
following observation, as detailed in the recent article by Bloch et. al [24]. They 
comment on the pressure-energy relation P = 2E/3 obeyed by a gas at unitarity. From 
this, it may be deduced that the effective two-body interaction potential has to be of 
inverse-square nature, given that it is not noninteracting. This also gives rise to the virial 
theorem given by Eq. (134) of their paper, which is in agreement with the experimental 
results of Thomas et. al. 1251. 
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